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Abstract 

In these four lectures I review the theory and phenomenology of weak decays 
of quarks, and their role in the determination of the parameters of the Stan- 
dard Model of particle physics, in testing subtle features of the theory and in 
searching for signatures of new physics. Attempts to understand CP- violation 
in current and future experiments is discussed. 

1 Lecture 1: Introduction 

Flavourdynamics, the study off the electroweak Lagrangian and its implications, is one of the central 
areas of research in particle physics. For example, in the United Kingdom the Mission Statement of the 
Particle Physics Community, as defined by our research council, is to obtain insights into the following 
three fundamental questions: 

• the origin of mass, 

• the 3 generations of elementary particles and their weak asymmetries, 

and 

• the nature of dark matter. 

The second of these topics is flavourdynamics, which is the subject of this lecture course. I will review the 
theoretical and experimental work being done in an attempt to determine the parameters of the Standard 
Model (SM) accurately, to test subtle properties of the SM and to define and search for signatures of 
New Physics. Although this primarily involves weak interactions, the major theoretical difficulty in 
interpreting experimental data on weak hadronic processes is controlling the non-perturbative strong 
interaction effects necessarily present in these decays. For this reason any review of flavour physics must 
contain a discussion also of these QCD effects, and much of the presentation below concerns this subject. 

In the next few years much of the effort of the experimental community will be concerned with 
trying to gain an understanding of CP- violation. Indeed, one of the Sakharov conditions for generating 
a matter-antimatter asymmetry in the universe in the big-bang cosmology is the requirement for the 
existence of CP-violation. The existence of three generations of quarks and leptons, and hence of a 
complex phase in the mixing matrix, implies the presence of CP- violation at some level. Very little 
is known, however, about the value of this phase and hence of the magnitude of CP-violation induced 
by this mechanism, and also about other possible sources. It is hoped that the intensive studies, which 
are about to begin will significantly further our understanding and this will be one of the main topics of 
discussion towards the end of this lecture course. 

The program for the four lectures is as follows: 

1 : The first lecture will be an introduction to weak decays in the Standard Model. Among the topics 
which will be covered will be charged and neutral currents, the Cabibbo-Kobayashi-Maskawa 
matrix and unitarity triangles, parity and charge conjugation symmetries, effective Hamiltonians 
and Operator Product Expansions and the Heavy Quark Effective Theory. 
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Fig. 1: The basic vertices representing the interactions of the quarks with the gauge bosons. The labels i and j represent the 
flavour quantum number (i, j = u,d,c,s, t, b). 



2: In the second lecture I will review the determination of the V c b and V u h matrix elements from 
semileptonic decays of 5-mesons. Lattice computations, which provide the opportunity to eval- 
uate non-perturbative strong interaction effects in weak decays in general, are introduced in this 
lecture. 

3: The third lecture will contain a review of K°-K° and B°-B° mixing. 

4: Finally, in the last lecture I review CP-violation in i?-decays and the subject of inclusive B- 
decays. 

The theoretical framework introduced in these lectures will be applied to several important physi- 
cal processes. It will not be possible, however, to discuss the full range of interesting processes which are 
providing, or will provide, fundamental physical information. A much larger set of physical quantities 
is considered in detail in the beautiful review by Buras and Fleischer [{j]], to which I refer the student, 
and from which I will quote extensively. I also refer the reader to ref. [§] for a broad introduction to 
the Dynamics of the Standard Model and to references [||] and [Q] for modern reviews of B-decays and 
CP-violation respectively. I will assume a familiarity with elementary Quantum Field Theory (as dis- 
cussed for example by J. Petersen and V. Zakharov at this school [^]), but will try to avoid technical 
complexities, focussing instead on the underlying ideas. 



1.1 The Interactions of Quarks and Gauge Bosons 

The interactions between the quarks and gauge bosons in the Standard Model are illustrated in Fig. [jj 
In these lectures we will be particularly interested in the weak interactions. The vertex for the charged 
current interaction in which quark flavour i changes to j is depicted in Fig 0(a) and has the Feynman rule 

i^^jT^C 1 - 76) , (!) 

where gi is the coupling constant of the SU (2) l gauge group and Vy is the ij element of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix (Vji = V*a) Eq. (Q) illustrates the V— A (vector— axial- vector) 
structure of the charged-current interactions. 

At low energies, so that the momentum in the W -boson is much smaller than its mass Myy, the 
four-fermion interaction mediated by the W -boson can be approximated by the local Fermi /3-decay 
interaction with coupling Gf, where 

G F g\ 
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(see Fig. ||). 
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Fig. 2: Approximation of the W-exchange interaction, by the four-fermion current-current vertex. 



1.2 CKM - Matrix 

The charged-current interactions are of the form 



J, 



cc 



(u,C,t)L7 At VcKM 




(3) 



where (assuming the Standard Model with 3 generations) the CKM-matrix, Vckm is a unitary 3x3 one, 
relating the weak and mass eigenstates. The 1996 particle data book [f7|] gives the following values for 
the magnitudes of each of the elements: 



/ 0.9745 - 0.9757 0.219 
0.218 - 0.224 0.9736 
v 0.004 - 0.014 0.034 



0.224 0.002 - 0.05 
0.9750 0.036 -0.046 
0.046 0.9989 -0.9993 



(4) 



The subscript L in eq. (|3|) represents left handed, (V>l = ^(1 — 7 5 )V0- 

If we have 2N g quark flavours, (N g is the number of generations) then Vckm is a N g x N g unitary 
matrix. It therefore has N 2 real parameters, however (2N g — 1) of these can be absorbed into unphysical 
phases of the quark fields FL leaving (N g — l) 2 physical parameters to be determined. 

In the two flavour case there is just one parameters, which is conventionally chosen to be the 
Cabibbo angle: 

/ cos C sin C \ 
\ — sin r cos r I 



(5) 



With three flavours there are 4 real parameters. Three of these can be interpreted as angles of ro- 
tation in three dimensions (e.g. the three Euler angles) and the fourth is a phase. The general parametri- 
sation recommended by the Particle Data Group |7|] is 



' C12C13 

-S12C23 - C12S23S13 exp(i<5i 3 ) C12C23 - 
\ S12S23 ~ ci 2 C23Si 3 exp(i(5i3) -C12S23 



S12C13 si 3 exp(-i(5i3j 
Sl2S23Sl3exp(i5i 3 ) S23C13 
" Sl2C 2 3Sl 3 exp(i5i 3 ) C23C13 



(6) 



where Cij and Sij represent the cosines and sines respectively of the three angles 6ij, ij = 12, 13 23. <5i 3 
is the phase parameter. 

The parametrisation in eq. @ is general, but awkward to use. For most practical purposes it 
is sufficient to exploit the empirical fact that the elements get smaller as one moves away from the 
diagonal of the matrix (see eq. (|j) ), and to use a simpler, but approximate parametrisation. It has become 
conventional to use the Wolfenstein parametrisation [M] : 
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'The reason why the number of phases is (2N g — 1) rather than 2N g , is that if the fields are multiplied by the same phase 
factor then J^ c is unchanged. Thus there is one fewer phase, which can be absorbed. 



A is approximately the Cabibbo angle, and the terms which are neglected in the Wolfenstein parametrisa- 
tion are of 0(A 4 ). Much of the current effort of the particle physics community is devoted to determining 
the four parameters A, A, p and rj, with ever increasing precision, and this will continue for some time to 
come. For this reason I will spend a significant fraction of these lectures on the theoretical issues related 
to the determination of these parameters. 

From the unitarity of the CKM-matrix we have a set of relations between the entries. A particularly 
useful one is: 

V ud V: b + V cd V: b + V td V t * b = . (8) 
In terms of the Wolfenstein parameters, the components on the left-hand side of eq. (S) are given by: 

AX 3 [p + if)} +0(A 7 ) 

-AX 3 + 0(A 7 ) (9) 
AX 3 [l-(p + i fj)] + 0(X 7 ), 

where p = p(l — A 2 /2) and f) = r/(l — A 2 /2). The unitarity relation in eq. (Js|> can be represented 
schematically by the famous "unitarity triangle" of Fig. |3| (obtained after scaling out a factor of ^4A 3 ). 
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Fig. 3: Unitarity Triangle corresponding to the relation in eq. (B). 



The notation in Fig. || is standard, and in the following lectures I will describe current and future 
attempts to determine the position of the vertex A, and hence the values of the angles a, (3 and 7. By 
using many different processes to overdetermine the position of A, one will be able to test the consistency 
of the SM, and check for evidence of physics beyond this model. 

1.2.1 The Cabibbo Sector 

Although much of the remainder of these lectures will be devoted to the determination of the entries 
in the CKM-matrix, since the Cabibbo sector is the best determined, I will not consider it beyond this 
subsection. I briefly review the status of the four entries in the 2x2 top-left submatrix of Vckm [0]- 

• V ud : 

I V ud \ = 0.9736 ±0.0010 . (10) 

This is the best determined of the elements and is obtained by studying superallowed /3-decays in 
nuclei (see fig. |](a)). Since the publication of the 1996 particle data book, a result with a smaller 
error has been presented, \V ud \ = 0.9740 ± 0.0005 [gj. 

• V us : 

\V US \ = X = 0.2205 ±0.0018 . (11) 




(a) (b) 




This matrix element is obtained from semileptonic decays of if -mesons (K + — ► 7r°e + ^ e , K® — > 
-K~e + v e , see fig. ||(b) ) or of hyperons. In the i'T-meson decays, parity symmetry implies that only 
the vector component of the V— A current contributes to the decays. Since the momentum transfer 
to the 7r-meson is small, and the vector current for degenerate quarks is conserved, the uncertainties 
due to strong interaction effects are small and can be estimated using chiral perturbation theory. 

V cd : 

\V cd \ = 0.224 ±0.016 . (12) 

The V c d element is obtained from charm production in deep inelastic neutrino (or antineutrino) 
nucleon scattering, see fig. ||(c). 

V cs : 

\Vca\ = 1.01 ±0.18 . (13) 

This matrix element can be determined from semileptonic decays of charmed mesons (see fig. 0(d)). 
In this case the large difference in the masses of the quarks make it difficult to estimate the strong 



interaction effects accurately (see sec. |2.1| ), and this is the reason for the relatively large error. 



The errors in eqs. (|l0[)-(|l"3"|) are those in the measurements of these matrix elements themselves. Com- 
paring these uncertainties with those in eq. (Q), where the constraints of unitarity have been imposed 
(assuming just three generations), we see how much these constraints reduce the errors. 



1.3 Flavour Changing Neutral Currents 

In the Standard Model, unitarity implies that there are no Flavour Changing Neutral Current (FCNC) 
reactions at tree level, i.e. that there are no vertices of the type: 

b s u c 



Quantum loops can, however, generate FCNC reactions, through box diagrams (see fig. |5]) or penguin 
diagrams (see fig. ^|), and we will discuss some of the physical processes induced by these loop-effects 
in the following lectures. 
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Fig. 5: Sox diagrams which contribute to the process of B°-B° mixing. 
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Fig. 6: Examples of penguin diagrams which contribute to the FCNC process d — » s. Diagrams (a) and (b) are electroweak 
penguins graphs, whereas diagram (c) is a gluonic penguin graph. 



As a consequence of the "GEVT'-mechanism [10] all the FCNC vertices vanish in the limit of 
degenerate quark masses. This is a consequence of the unitarity of the CKM-matrix. For example, 
consider the vertex in the diagram of fig. |6|(a), in the hypothetical situation in which there is a "horizontal" 
symmetry such that m u = m c = m t . In this case the contribution from each of the positively charged 
quarks in the loop is simply proportional to the corresponding CKM-matrix elements, so that the total 
contribution is 

V ud V* s + VcdV*s + V t dV t * s (14) 
which vanishes by unitarity of the CKM-matrix. A similar argument holds for all the FCNC vertices. 

Of course, in reality the masses of the quarks are not equal and the FCNC vertices are not zero. 
Depending on the process being studied, however, the GIM-mechanism may lead to a suppression of the 
amplitude. 

1.4 P, C and CP Symmetries 

Symmetries play a fundamental role in physics in general and in particle physics in particular. In this 
subsection I briefly introduce the violation of CP-symmetry, which will be a focus of many of the studies 
at P-physics facilities in the next few years. 

Parity: The parity transformation reverses the sign of spatial coordinates, (x, t) — ► (— x, t). The vector 
and axial- vector fields transform as: 

VJx,t) -> V^(-x,t) and AJx,t) -A"(-x,t) , (15) 



and the vector and axial- vector currents transform similarly. Left-handed components of fermions (i (1 — 
7 5 ) , 0) transform into right-handed ones (|(1 + 7 5 )V0> and vice-versa. Since weak interactions only 
involve the left-handed components, parity is not a good symmetry of the weak force, in distinction to 
the strong and electromagnetic forces (QCD and QED are invariant under parity transformations). 



Charge Conjugation: For free fields (4>(x) say) we can perform a fourier decomposition, with coeffi- 
cients which contain annihilation and creation operators: 

00) = ••• a + ■■■ , (16) 

where a and represent the annihilation and creation operators for the particle and antiparticle. The 
charge conjugation transformation is the interchange of a and b. In an interacting field theory this is not 
in general the same as interchanging physical particles and antiparticles for which we need the CPT 
combined transformation (T is time reversal transformation). Nevertheless C-tranformations are also 
interesting in field theories. Under C the currents transform as follows: 

^17/^2 -> -4>2l^i and ^17^75^2 -> ^27^75^1 , (17) 
where represents a spinor field of type (flavour or lepton species) i. 



CP : Under the combined CP-transformation, the currents transforms as: 

ipi7^ip2 -> -^27^1 and ^17^75^2 -> -^27^75^1 , (18) 

where the fields on the left (right) hand side are evaluated at (x, t) ( (— x, t) ). Consider now a charged 
current interaction: 

(K ~ iW D ^V(l - l b )VijD j + {Wl + iWl) 5V(1 - 7 5 )^ *U { , (19) 
where W represents the field of the VF-boson and U l and are up and down type quarks of flavours i 



and j respectively. Under a CP transformation, the interaction term in eq. ( |19| ) transforms to: 

(K + iwi) - ^VijU' + (wi - iwi) wra - i 5 m*D j , (20 



where the parity transformation on the coordinates is implied. Comparing eqs. ( |l9| ) and (20), we see 
that CP-conservation requires the CKM-matrix Vij to be real (or more strictly that any phases must be 
absorbable into the definition of the quark fields). If the only source of CP-violation in nature is the 



phase in the CKM-matrix, then at least three generations of quarks are required (see subsec. L2 ). A 
central element of the research of the forthcoming generation of experiments in P-physics will be to 
determine whether the phase in the CKM-matrix is the only (or the main) source of CP-violation. 

CP- violation appears to be small compared to the strength of the weak interaction, so that CP is 
a fairly good symmetry. We recall that the presence of CP- violation is one of the Sakharov conditions 
for the creation of the baryon-antibaryon asymmetry of the universe. 



1.5 Leptonic Decays of Mesons 

As already mentioned, the main difficulty in making predictions for weak decays of hadrons is in con- 
trolling the non-perturbative strong-interaction effects. We will encounter this problem several times in 
subsequent lectures, but, as an introduction, I now discuss it briefly in a particularly simple situation, that 
of the leptonic decays of pseudoscalar mesons in general, and of the P-meson in particular. 
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Fig. 7: Diagram representing the leptonic decay of the B-meson. 



The leptonic decay of a pseudoscalar meson (e.g. the i?-meson) is illustrated in the diagram of 
Fig. 0. All the QCD effects corresponding to the strong interactions between the quarks in Fig. [7] are 
contained in the matrix element of the V— A hadronic current: 

(0\b^(l-l 5 )u\B(p)} . (21) 

Lorentz and parity symmetries imply that the matrix element of the vector component of the V— A 
current vanishes ((0| b^u \B(p)) = 0). This component is an axial-vector, and there is no axial-vector 
which can be constructed from the single vector we have at our disposal (p) and the invariant tensors 
(g IM ' u and e^ upu ). These symmetries also imply that the matrix element of the axial component of the 
V— A current is a Lorentz- vector, and hence that it is proportional to 

(0\bj^ 5 u\B(p)} =if B p fi - (22) 

The constant of proportionality, /g, is called the decay constant of the 5-meson. Thus all the QCD- 
effects for leptonic decays of B-mesons are contained in a single number, fg. An identical discussion 
holds for other pseudoscalar mesons (tt, K, ■ • •)■ The normalisation used throughout these lectures, corre- 
sponds to f w ~ 132 MeV (/b is unknown at present). For leptonic decays of vector mesons, space-time 
symmetries also imply that the strong interaction effects are contained in a single decay constant. Thus, 



the role of quantitative non-perturbative methods, such as lattice QCD (discussed in sec. 2.3 below) or 



QCD sum-rules, is to provide the framework for the evaluation of matrix elements such as that in eq. (pl|). 
1.6 Operator Product Expansions and Effective Hamiltonians 

Even though this is a course in flavour physics, we cannot avoid the fact that the quarks are interacting 
strongly and hence that QCD effects must be considered. Indeed, as already mentioned above, our 
rather primitive control of non-perturbative QCD is the principal source of uncertainty in interpeting 
experimental data of weak decays and in determining the fundamental parameters of the Standard Model. 

In Paolo Nason's lecture course [p|], we have seen how the property of asymptotic freedom, which 
states that the interactions of quarks and gluons become weaker at short distances, enables us to use 
perturbation theory to make predictions for a wide variety of short-distance (or light-cone) dominated 
processes. For separations \x\ <C ^-qcd (M < 0-1 f m sav ) or corresponding momenta \p\ 3> Aqcd 
(\p\ > 2 GeV say), we can use perturbation theory. The natural scale of strong interaction physics is 
of 0(1 fm), however, and so in general, and for most of the processes discussed in this course, non- 
perturbative techniques must be used. 

As an example, consider a decay of a ET-meson into two pions, for which a tree-level diagram of 
the quark subprocess is depicted in Fig. ||(a). The amplitude is proportional to 

^| KdYus (7T7r|(dy*(l - j 5 )u) (07^(1 - 7 5 )s)\K) , (23) 

i.e. it is proportional to the matrix element of the operator 



Oi = (J 7 M (1 - 7 5 K) (n 7M (l-7 5 )s) 



(24) 
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Fig. 8: Diagrams contributing to hadronic decays of AT -mesons 



between \K) and \tt-k) states. A one-loop (i.e. one-gluon exchange) correction to this diagram is shown 
in Fig. |8|(b). This clearly generates a second operator, (jT a 7^(l — 7 5 )it) (uT a ^f fM (l — 7 5 )s) [], which, 
using Fierz identities can be written as a linear combination of 0\ and O2 where 

2 = (d>f(l - 7 5 )s) (« 7m (1 - 7 5 )«) • (25) 

This discuusion can be generalized to higher orders, as I now explain. 

Using the operator product expansion (OPE) one writes the amplitude for a weak decay process 
from an initial state X to a final state T in the form: 

T TI = % Vckm £ a (fx) (T\ Oi Ou) I J) . (26) 
v2 i 

fi is the renormalization scale at which the composite operators Oi are defined and Vck M lS the appro- 



priate product of CKM-matrix elements. The expansion (g6j) is very convenient. The non-perturbative 
QCD effects are contained in the matrix elements of the operators Oi, which are independent of the large 
momentum scale, in this case of M\y- The Wilson coefficient functions Ci(/j,) are independent of the 
states Z and T and can be calculated in perturbation theory. Since physical amplitudes manifestly do not 
depend on fi, the ^-dependence in the operators Oi(fi) is cancelled by that in the coefficient functions 
Ci{n). The effective Hamiltonian for weak decays can then be written in the form: 

H cS = % Vckm £ d ( M ) 0< (/*) . (27) 
V2 



In order to gain a little further intuition into the origin of eq.(|26|) consider the one-loop graph of 
Fig. ||(b). Dimensional counting at large loop momenta k readily shows that the graph is ultra-violet 
convergent (in the Feynman gauge, say): 

k large k k k 2 k 2 — Myy ' v ' 

where the first two factors in the integrand (1/k) correspond to the two quark propagators, the third 
(l/k 2 ) to the gluon propagator and the final factor to the propagator of the W -boson. It can be deduced 
from eq.(|2^) that the contribution from this graph contains a term proportional to log(M^/p 2 ), where p 
is some low-momentum (infra-red) scale. 

Consider now the one-loop graph of Fig. ^, depicting a correction to a matrix element of the 
operator 0\ defined in eq. ([24)) . Now the VF-propagator is absent, and the power counting in the large 
momentum region gives: 

11 L d*k, (29) 



I k large k k k 2 

2 The T a 's are the eight matrices, representing the generators of the SU (3) colour group in the fundamental representation. 
The coupling of quarks and gluons is proportional to the elements of these matrices. 



Fig. 9: One-loop correction to the Green function of the operator 0\. 



i.e. a logarithmically divergent contribution. After renormalization at a scale \i there will be a contribu- 
tion to this graph proportional to log(// 2 /p 2 ). In the low momentum (non-perturbative) region (i.e. for 
\hp\ *C M\y) the contributions from the graphs in Figs. |8|(b) and |9] are the same. Thus we can split the 
factor of log(Mjy /p 2 ) from the graph of Fig. [|(b) into two parts; the first, log(fi 2 /p 2 ) is contained in the 
matrix element of 0\, and the second, log(M^/// 2 ), into the coeficient function C\. Since the infra-red 
behaviour of the graphs in Figs, ^b) and |9] are the same, the coefficient functions are independent of the 
treatment of the low-momentum region. Moreover, in order to calculate the coefficient functions, we can 
choose any convenient external states, and in practice one chooses quark or gluon states. We refer to the 
evaluation of the coefficient functions as the process of matching. 

In iT-decays it is natural, although not necessary, to choose fj, to be as small as one can without 
invalidating perturbation theory. In this way one might hope to avoid large logarithms (log(/i 2 /p 2 ) ) 
which would spoil any insights we might have about the operator matrix elements from non-relativistic 
quark models or other approaches. Of course this implies the presence of large logarithms of the type 
a,g \og n {My V / fj, 2 ) in the coefficient functions, but these can be summed using the renormalization group, 
leading to factors of the type: 

-a s (M w )^°W° 



, ] (30) 

where 70 is the one-loop contribution to the anomalous dimension of the operator (proportional to the 
coefficient of log(/j, 2 /p 2 ) in the evaluation of the graph of fig. |8^b) ) and j3o is the first term in the (3- 
function, ((3 = dg/ 91n(//) = — (3q g 3 /\67r 2 ) []. Eq. (3C) represents the sum of the leading logarithms, 



i.e. the sum of the terms a™ log n (M^//i 2 ). For almost all the important processes, the first (or even 
higher) corrections have also been evaluated. 

We end this subsection by making some further points about the use of the OPE in weak decays 
and other hard processes: 

• We shall see below that for some important physical quantities (e.g. e'/e), there may be as many 
as ten operators, whose matrix elements have to be estimated. 

• One may try to evaluate the matrix elements using some non-perturbative computational technique, 
such as lattice simulations or QCD sum-rules (see below), or to determine them from experimental 
data. In the latter case, if there are more measurements possible than unknown parameters (i.e. 
matrix elements), then one is able to make predictions. 

• In weak decays the large scale, My/, is of course fixed. For other processes, most notably for deep- 
inelastic lepton-hadron scattering, the OPE is useful in computing the behaviour of the amplitudes 
with the large scale (e.g. with the momentum transfer). 



3 In general when there is more than one operator contributing to the right hand side of eq. (^), the expression in eq. (|3C 
generalized into a matrix equation, representing the mixing of the operators 



1.7 The Heavy Quark Effective Theory - HQET 



Most of the important physical properties of the hydrogen atom are independent of the mass of the 
nucleus. Analogous features hold in the physics of heavy hadrons, where by heavy I mean that mq 3> 
Aqcd and tuq is the mass of the heavy quark. During the last few years the construction and use of the 
Heavy Quark Effective Theory (HQET), has proved invaluable in the study of heavy quark physics. In 
this subsection I briefly introduce the main features of the HQET, which will be used in the following 
lectures. 

Consider the propagator of a (free) heavy quark: 

If the momentum of the quark p is not far from its mass shell, 

Pn = rriQVn + , (31) 



(v 2 = 1), then 



where \k^\ -C tuq and is the (relativistic) four velocity of the hadron containing the heavy quark 

► = j kkJl i + p(\M\ 

p 2 v-k+ie \ tuq ) " 

(1+ y)/2 is a projection operator, projecting out the large components of the spinors. This propagator 
can be obtained from the gauge-invariant action 

Chqet = h(iv • D)^-h , (32) 

where h is the spinor field of the heavy quark Q Chqet is independent of tuq, which implies the exis- 
tence of symmetries relating physical quantities corresponding to different heavy quarks (in practice the 
b and c quarks). The light degrees of freedom are also not sensitive to the spin of the heavy quark, which 
leads to a spin-symmetry relating physical properties of heavy hadrons of different spins. Consider, for 
example, the correlation function: 

J d 3 x(0\J H (x)J H (0)\0) , (33) 

where J H (Jh) is an interpolating operator which can create a heavy hadron H, which we take to be a 
pseudoscalar or vector meson. The hadron is produced at rest, with four velocity v = (0, 1). We can use 
the interpolating operator Jh = hj 5 q for the pseudoscalar meson and Jh = hrfq (i = 1, 2, 3) for the 
vector meson. This means that the correlation function will be identical in the two cases except for the 
factor 

.1+7° K 1-7° 
7 5 ^7 5 = —2- (34) 

when if is a pseudoscalar meson, and 

■ 1 + 7° ■ 1-7° 

y^y = -3^ (35) 

when it is a vector meson. Since the correlation functions behave with time as exp(— iMfft) , this implies 
that the pseudoscalar and vector mesons are degenerate, up to relative corrections of 0(AQ CD /m,Q): 

M P = M V + 0(A 2 QCD /m Q ) , (36) 



4 Of course there are formal derivations of the action in eq. 



where Mp and My are the masses of the pseudoscalar and vector mesons respectively, or equivalently 

My - Mp = constant . (37) 



This relation is reasonably well satisfied experimentally: 

Af|. - M% = 0.485 GeV 2 and M 2 D , -M 2 D = 0.546 GeV 2 , (38) 

which is encouraging. However, I also mention in passing that the light mesons, to which the HQET 
certainly does not apply, also satisfy these relations numerically: 

M\* — M\ = 0.552 GeV 2 and M 2 - M 2 = 0.571 GeV 2 . (39) 

The usefulness of the HQET does not lie in the fact that the residual momenta k are always much 
smaller than Aqcd', indeed this is not true, in QCD one does have hard gluons with arbitrarily large mo- 
menta. However the effects of hard gluons can be evaluated in perturbation theory. The non-perturbative 
strong interaction effects are the same in QCD and in the HQET (up to corrections of 0(Aqcd A"q)> 
which for the moment we neglect), and so the HQET relations, such as that between the decay constants 
of pseudoscalar and vector mesons, which can be deduced from the above proportionality of correlation 
functions, are violated only by perturbatively calculable corrections. 

The principal application of the HQET is in 6-physics, and to a lesser extent to charm physics. 
In both cases the corrections of 0(Aqc£)/?tiq) are significant, and hence one would like to compute 
these corrections. In practice, this gets progressively more difficult, and it is a much debated point as to 
whether even the first corrections (i.e. those of 0{KQCD/fnQ) relative to the leading terms) have been 
reliably calculated for any interesting quantity. 



2 Lecture 2: V c b and V u b 

In this lecture I will discuss exclusive semileptonic decays of S-mesons, in which the fc-quark decays 
into a c- or u-quark, and from which one can determine the V c b and V u b elements of the CKM-matrix. 
The decays are represented in Fig. |K]. It is convenient to use space-time symmetries to express the matrix 
elements in terms of invariant form factors (using the helicity basis for these as defined below). When 
the final state is a pseudoscalar meson P, parity implies that only the vector component of the V— A 
weak current contributes to the decay, and there are two independent form factors, / + and f°, defined 
by 



(P(k)\V fl \B(p))=f + (q 2 



m r> — mi 



+/V) B (40) 



q 2 



where q is the momentum transfer, q = p—k. When the final-state hadron is a vector meson V, there are 
four independent form factors, V, A\, A2 and A: 

(V(k,e)\V»\B(p)) = 2V{q2) e^efakg (41) 
mp + ray M 

(V{k,e)\A"\B{p)) = i(m B+mv )A 1 {q 2 )e^ -i A ^ e*-p{p+kf 

mp+mv 

A(a 2 ) 

+ i^-^2myeW , (42) 

where e is the polarization vector of the final-state meson, and q = p—k. Below we shall also discuss 
the form factor Aq, which is given in terms of those defined above by Aq = A + A3, with 

A 3 = m B + my M _ mp-my^ 
2mv 2my 



V—A 




D, D*, TV, p 



Fig. 10: Diagram representing the semileptonic decay of the B-meson. q represents the light valence antiquark, and the black 
circle represents the insertion of the V—A current with the appropriate flavour quantum numbers. 



2.1 Semileptonic B —> D and B —> D* Decays 

Semileptonic B — > D* and, more recently, B — > D decays are used to determine the V c b element of the 
CKM matrix. Heavy quark symmetry is rather powerful in controlling the theoretical description of these 
heavy-to-heavy quark transitions (see e.g. the review article by Neubert for details and references [12]). 

In the heavy quark limit, all six form factors in eqs. ( |40| ) - ( f42| ) are related and there is just one uni- 
versal form factor known as the Isgur-Wise (IW) function, which contains all the non-perturbative 
QCD effects. Specifically: 



/V) = V(q 2 ) = A (q 2 ) = A 2 (q 2 ) 



(m B + m D ) 2 



A,{q 2 



ms + mo 
2^niBm D 



£M, (44) 



where u> = vb-vd is the velocity transfer variable (vb,d are the four velocities of the corresponding 
mesons). Here the label D represents the D- or D*-meson as appropriate (pseudoscalar and vector 
mesons are degenerate in this leading approximation). Vector current conservation implies that the IW- 
function is normalized at zero recoil, i.e. that £(1) = 1. This property is particularly important in the 
extraction of 

To get some insights into the origin of the heavy quark relations above, let us consider B — > D 
decays. We can rewrite eq. (4C) in terms of the form factors / + and /~ (instead of using the helicity 
basis as in eq. (HO) ): 



(D(p')\c^b\B( P )) = f + (q 2 )(p+ P r - r(q 2 )(p 



p 



(45) 



where q = p — p'. Defining the four- velocities v and v' by p = tub v and p' = vtld v', eq. (45) can be 
rewritten as: 



(D(v')\c^b\B(v)} 



(m B + m D )f + {q 2 ) - (m B - m D )f (q 2 )] (v + v'Y 



+ - 



(m B - m D )f + (q 2 ) - (m B + m D )r(q 2 )} (v - v'Y • (46) 



Now in the HQET: 



^m B m D (P(v')\ Kn^K \P(v)) = ^/m B m D (,(uj) (v + v'Y , 



(47) 



where uj = v-v', and the projection operators (e.g. (l+^)/2) have been absorbed into the heavy 
quark fields. P represents a heavy pseudoscalar meson, and the factors ^/mBrnB are chosen in or- 
der to introduce a mass-independent normalization of states (so that the matrix elements are independent 
of any masses). There is no term proportional to (v 
h v i(tf— rf')h v = 0, and so we obtain the relations: 



on the right hand side of eq. (^), since 



r 



m B - m D 
iris + mn 



f+ and / 



± _ m B ± m D 



CM 



(48) 



To understand the normalization condition £(1) = 1, consider the forward matrix element: 



(B(p)\ b^b \B(p)) = 2pff + (0) = 2pf , (49) 

where the last relation follows from the conservation of the current b^b which implies that / + (0) = 1. 
In the HQET this matrix element is equal to tub £(1), and so £(1) = 1. The spin symmetry relates 
the form-factors of B — > D* decays to 

The relations in eq. ( (44| ) are valid up to perturbative and power corrections. The precision with 
which can be extracted is limited by the theoretical uncertainties in the evaluation of these corrections. 
Nevertheless we are in the fortunate situation that it is uncertainties in corrections (which are therefore 
small) which control the precision. 

Allowing for these corrections to the results in the heavy quark limit, one writes the decay distri- 
bution for B —y D* decays as: 



dT G 



2 



du 48?r 3 



^(m B -m D *) 2 m 3 D * \Aj 2 -1 (w+1) 2 



^ 4a; m B — 2u:mBrnr)* + m 2 D * 



\V cb \ 2 T 2 {u) 



uj+1 (niB — rnr)*) 2 

(50) 

where is the IW-f unction combined with perturbative and power corrections. It is convenient 

theoretically to consider this distribution near u> = 1, in order to exploit the normalization condition 
£(1) = 1. In this case there are no 0(1/toq) corrections (where Q = b or c) by virtue of Luke's 



theorem []13[], so that the expansion of ^(l) begins like: 



m=riA[l + 0^+c 2 ^ + ...) , (51) 
rriQ m£ 



where t)a represents the perturbative corrections. The one-loop contribution to t\a has been known for 
some time now, whereas the two-loop contribution was evaluated last year, with the result [|l4|] : 

t)a = 0.960 ± 0.007 , (52) 

where I have taken the value of the two loop contribution as an estimate of the error. 

The power corrections are much more difficult to estimate reliably. Neubert has recently combined 
the results of refs. [|i"5|, [l6|, [i7|] to estimate that the 0(l/mg) terms in the parentheses in eq. ( |5l| ) are about 
-0.055 ± 0.025 and that = 0.91(3) . Bigi, Shifman and Uraltsev 0, consider the uncertainties 
to be bigger and obtain 0.91(6). Combining the latter, more cautious, theoretical value of .F(l), with the 



experimental result [|9|] .F(l)|V c6 | = (34.3±1.6)10~ 3 , readily gives |V c6 | = (37.7±1.8 exp ±2.5 th )10- 3 . 



J exp- 

Here I am only discussing exclusive decays. Buras and Fleischer, perform an analysis of both 
exclusive and inclusive semileptonic decays and quote V c b = (40 ± 3) 10~ 3 as their best value [jl|]. 

IV^I gives the scale of the unitarity triangle: 

X\V cb \ = X 3 A. (53) 



2.2 Semileptonic B —> ir and B — > p Decays 

In this subsection we consider the heavy-to-light semileptonic decays B — > p and B —> it which are now 
being used experimentally to determine the V u b matrix element [2C, 21). Heavy quark symmetry is less 
predictive for heavy-to-light decays than for heavy-to-heavy ones. In particular, there is no normalization 
condition at zero recoil corresponding to the condition £(1) = 1, which is so useful in the extraction of 
V c b (see subsection |Q| ). The lack of such a condition puts a premium on the results from nonperturbative 
calculational techniques, such as lattice QCD or light-cone sum rules. Heavy quark symmetry does, 



however, give useful scaling laws for the behaviour of the form factors with the mass of the heavy quark 
(rriQ) at fixed u: 



f + ,A ,A 2 ,V~ 



1 



^3 ~ m j 



(54) 



These scaling relations are particularly useful in lattice simulations, where the masses of the quarks are 
varied. Moreover, the heavy quark spin symmetry relates the B — > V matrix elements [22, |3|] (where 



V is a light vector particle) of the weak current and magnetic moment operators, thereby relating the 
amplitudes for the two physical processes B° — ► p + l~i>i and B — ► K*j, up to SU (3) flavour symmetry 
breaking effects. These relations also provide important checks on theoretical, and in particular on lattice, 
calculations. 

Recent compilations of results for iKb/Kftl, include 0.06-0.11 obtained from exclusive decays 
only (or 0.08 ± 0.01 eX p ± 0.02 th from inclusive decays) [Q and 0.08 ± 0.02 from both inclusive 



and exclusive decays [j^]. In terms of the parameters of the unitarity triangle (see fig. ^]): 



CA 



\Vu d v: b \ 



p2 + fj2 = (1 



2 ' 



1 \V ub \ 



AAA 



I Kb | 



(55) 



\VcdV^\ vr ' v 2 'X \V cb \ ' |y c6 | 

Thus a measurement of |Kb/^c&| implies that the locus of possible positions of the vertex A is a circle 
centred on C (see fig. |ll|). In practice, of course, because of the theoretical and experimental uncertain- 
ties, the circle is replaced by a band. 



A(p,fj) 




C(0,0) 



5(1,0) 



Fig. 11: A precise determination of |Vui,/\4i,| would fix the vertex A to lie on a given circle centered on C (schematically 
represented by the solid line). 



2.3 Lattice QCD 

We have seen that the main difficulty in making predictions for weak hadronic decays is our inability to 
control the long-distance strong interaction effects present in these decays. The principal technique for 
the evaluation of these non-perturbative effects is lattice QCD, and I now digress from the main theme 
of this course to outline briefly the principles underlying lattice calculations of non-perturbative QCD 
effects in weak decays and to discuss the main sources of uncertainty present in these computations. 

In subsection we have seen that, by using the Operator Product Expansion, it is generally 
possible to separate the short- and long-distance contributions to weak decay amplitudes into Wilson 
coefficient functions and operator matrix elements respectively. Thus, in order to evaluate the non- 
perturbative QCD effects, it is necessary to compute the matrix elements of local composite operators. 
This is achieved in lattice simulations, by the direct computation of correlation functions of multi-local 
operators composed of quark and gluon fields (in Euclidean space): 

(0|O(xi,x 2 ,...,^ n )|0) = ^J[DA fM \[Dip}[D$]e- s O(x 1 ,x 2 ,...,x n ), (56) 
where Z is the partition function 

Z = JpA^LhfllD-fie- 3 , (57) 



S is the action and the integrals are over quark and gluon fields at each space-time point. In eq. (56) 
0(xi,X2, ■ ■ ■ , x n ) is a multi-local operator; the choice of O governs the physics which can be studied. 

We now consider the two most frequently encountered cases, for which n=2 or 3. As a first 
example, let 0(x±, x 2 ) be the bilocal operator 

2 (x 1 ,x 2 ) = T{J H (x 1 )j\ I (x 2 )} , (58) 

where Jh is an interpolating operator for the hadron H whose properties we wish to study. In lattice 
computations we evaluate the two point correlation function 

C 2 {t x ) = ^e i P- x (0|O 2 (x,0)|0) . (59) 

X 

Inserting a complete set of states {\n}} we have: 

C2(t x ) = EE e ip - x (0\J H (x,t)\n) (n | 4(0,0) |0> , (60) 

n x 

= J2 etP ' X (°\ J H(^t)\H)(H\ 4(0,0)|0) +•••. (61) 

X 

In eq. (|6l]), the ellipsis represents the contributions from heavier states than H, which we assume to be 
the lightest hadron which can be 
argument of Jh to zero, we find I 



the lightest hadron which can be created by the operator 4- Using the translation operator to move the 



C2(t x ) = ^ w -\(0\J H (0)\h)\ 2 + •■• , (62) 

where Eh = J M H + p 2 . At large positive times t x the contribution from each heavier hadron, H' with 
mass mu' say, is suppressed by an exponential factor, exp ( — (Eh> — En)t x ), so that the contribution 
from the lightest state H is isolated. This is illustrated in the following diagram: 

H 

t x ■ 

In lattice simulations the correlation function C 2 is computed numerically, by discretizing space- 
time (hence the word lattice), evaluating the functional integral in eq. (|56|) by Monte-Carlo integration. 
By fitting the results to the expression in eq. (62) both the mass Mh and the matrix element ( | Jh (0) | h ) 
can be determined []. 

As an example consider the case in which H is the S-meson and Jh is the axial current (with 
the flavour quantum numbers of the S-meson). In this case one obtains the value of the leptonic decay 
constant Jb, defined in eq. (f?2|). 

It will also be useful to consider three-point correlation functions: 

x,y 

where, J\ and J 2 are the interpolating operators for hadrons H\ and H 2 respectively, O is a local operator, 
and we have assumed that t x > t y > 0. Inserting complete sets of states between the operators in eq. ( |63| ) 
we obtain 

e -E!ty e -E 2 (t x -ty) 

1E\ 2E 2 

(H 2 (p, E 2 )\O(0, 0)|ffi (p+q, Ei)) (ffi(p+q, 4(0, 0)|0) + • • • , (64) 

s The phase-space factor of \I2Eh in eq. ( |62| ) is implicitly included as part of the definition of the summation over states in 
eqs. @ and @. 

frequently it is most convenient to evaluate the correlation function with the hadron at rest, i.e. with p = 0. 



C S (t x ,t y ) = E e ip - x e^ y (0 | J 2 (x,t x )0(y,t y ) 4(0,0)10) , (63) 



C 3 (t x ,t y ) = — — (0|J 2 (0,0)| J f7 2 (p, J B 2 ))x 



where E\ = y Mf Ii + (p+q) 2 , Ei = yM| 2 + p 2 and the ellipsis represents the contributions from 
heavier states. The exponential factors, exp(— E\t y ) and exp ( — E2(t x — t y )), ensure that for large 
time separations, t y and t x — t y , the contributions from the lightest states dominate. The three-point 
conelation function is illustrated in the diagram 

H x O H 2 
• ► • ► • 

t y t x 

All the elements on the right-hand side of eq.(|54|) can be determined from two-point correlation func- 
tions, with the exception of the matrix element (H2\0\H\). Thus by computing two- and three-point 
correlation functions the matrix element (H2\0\H\) can be determined. 

The computation of three-point correlation functions is useful, for example, in studying semilep- 
tonic and radiative weak decays of hadrons, e.g. if H\ is a i?-meson, H2 a D meson and O the vector 
current Iry^c, then from this correlation function we obtain the form factors relevant for semileptonic 
B —> D decays. 

I end this brief summary of lattice computations of hadronic matrix elements with a word about the 
determination of the lattice spacing a. It is conventional to introduce the parameter (3 = Q/g^ (a) , where 
go (a) is the bare coupling constant in the theory with the lattice regularization. It is j3 (or equivalently 
go(a)) which is the input parameter in the simulation, and the corresponding lattice spacing is then 
determined by requiring that some physical quantity (which is computed in lattice units) is equal to the 
experimental value Q For example, one may compute m p a, where m p is the mass of the p-meson, and 
determine the lattice spacing a by dividing the result by 769 MeV. 



L 

Fig. 12: Schematic diagram representing a lattice containing a hadron. a and L are the lattice spacing and length of the lattice 
respectively. 



2.3.1 Sources of Uncertainty in Lattice Computations: 

Although lattice computations provide the opportunity, in principle, to evaluate the non-perturbative 
QCD effects in weak decays of heavy quarks from first principles and with no model assumptions or free 
parameters, in practice the precision of the results is limited by the available computing resources. For 
these computations to make sense it is necessary for the lattice to be sufficiently large to accommodate 
the particle(s) being studied (L » 1 fm say, where L is the spatial length of the lattice), and for the 
spacing between neighbouring points (a) to be sufficiently small so that the results are not sensitive to 

7 The bare quark masses are also parameters which have to be determined; one has to use as many phyical quantities as there 
are unknown paramters. 



the granularity of the lattice (oAqcd <C 1), see fig. 12. The number of lattice points in a simulation is 



limited by the available computing resources; current simulations are performed with about 16-20 points 
in each spatial direction (up to about 64 points if the effects of quark loops are neglected, i.e. in the so 
called "quenched" approximation). Thus it is possible to work on lattices which have a spatial extent of 
about 2 fm and a lattice spacing of 0. 1 fm, perhaps satisying the above requirements. I now outline the 
main sources of uncertainty in these computations: 



Statistical Errors: The functional integrals in Eq. ( p6[ ) are evaluated by Monte-Carlo techniques. 
This leads to sampling errors, which decrease as the number of field configurations included in the 
estimate of the integrals is increased. 

Discretization Errors: These are artefacts due to the finiteness of the lattice spacing. Much effort 
is being devoted to reducing these errors either by performing simulations at several values of the 
lattice spacing and extrapolating the results to the continuum limit (a = 0), or by "improving" 
the lattice formulation of QCD so that the error in a simulation at a fixed value of a is formally 



smaller [ 25 ] — [ 28 ] . In particular, it has recently been shown to be possible to formulate lattice QCD 



in such a way that the discretization errors vanish quadratically with the lattice spacing 0290, even 



for non-zero quark masses [T30J]. This is in distinction with the traditional Wilson formulation Q31[], 
in which the errors vanish only linearly. In some of the simulations performed in recent years 
improved actions and operators have been used. In some of these these studies the improvement 
is implemented at tree-level, so that the artefacts formally vanish more quickly (like aa s (a)) than 
for the Wilson action. This tree-level improved action is denoted as the SW (after Sheikholeslami- 



Wohlert who first proposed it [27]) or "clover" action 



Extrapolations to Physical Quark Masses: It is generally not possible to use physical quark masses 
in simulations. For the light (u- and d-) quarks the masses must be chosen such that the correspond- 
ing 7r-mesons are sufficiently heavy to be insensitive to the finite volume of the lattice. In addition, 
as the masses of the quarks are reduced the computer time necessary to maintain the required level 
of precision increases rapidly. For the heavy quarks Q (i.e. for c, and particularly for b) the masses 
must be sufficiently small so that the discretization errors, which are of 0(m,Qa) or 0(rriQa 2 ), are 
small. The results obtained in the simulations must then be extrapolated to those corresponding to 
physical quark masses. 

Finite Volume Effects: We require that the results we obtain be independent of the size of the lattice. 
Thus, in principle, the spatial size of the lattice L should be 3> 1 fm (in practice L > 2-3 fm), and 
as mentioned above, we cannot use very light u- and ci-quarks (in order to avoid very light pions 
whose correlation lengths, i.e. the inverses of their masses, would be of O(L) or greater). 

Contamination from Excited States: These are the uncertainties due to the effects of the excited 



states, represented by the ellipsis in Eq. (p2|). In most simulations, this is not a serious source 
of error (it is, however, more significant in computations with static quarks). Indeed, by evalu- 
ating correlation functions (Jjj(x) J' H (0)} for a variety of interpolating operators {Jh, J'h}> it is 
possible to obtain the masses and matrix elements of the excited hadrons (for a recent example 
see [|2|]). 

Lattice-Continuum Matching: The operators used in lattice simulations are bare operators de- 
fined with the lattice spacing as the ultra-violet cut-off. From the matrix elements of the bare 
lattice composite operators, we need to obtain those of the corresponding renormalized operators 
defined in some standard continuum renormalization scheme, such as the MS scheme. The re- 
lation between the matrix elements of lattice and continuum composite operators involves only 
short-distance physics, and hence can be obtained in perturbation theory. Experience has taught 
us, however, that the coefficients in lattice perturbation theory can be large, leading to significant 
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Fig. 13: Comparison of lattice spacings (the quantity plotted is actually the inverse lattice spacing, a' 1 ) determined from dif- 
ferent physical quantities in quenched and unquenched simulations by the SESAM and T^L collaborations F denotes the 
value determined from the static quark potential and p denotes the value determined from the p-meson mass, while 2S-1S and 
IP-IS denote values obtained from energy level splittings in QQ bound states using the lattice formulation of nonrelativistic 
QCD. 



uncertainties (frequently of 0(10%) or more). For this reason, non-perturbative techniques to eval- 
uate the renormalization constants which relate the bare lattice operators to the renormalized ones 



have been developed, using chiral Ward identities where possible [33] or by imposing an explicit 



renormalization condition [34] (see also Refs. [29, 35]), thus effectively removing this source of 
uncertainty in many important cases. 

"Quenching": In most current simulations the matrix elements are evaluated in the "quenched" 
approximation, in which the effects of virtual quark loops are neglected. For each gluon configu- 



ration {Ufj,(x)}, the functional integral over the quark fields in Eq. (56) can be performed formally, 
giving the determinant of the Dirac operator in the gluon background field corresponding to this 
configuration. The numerical evaluation of this determinant is possible, but is computationally 
very expensive, and for this reason the determinant is frequently set equal to its average value, 
which is equivalent to neglecting virtual quark loops. Gradually, however, unquenched calcula- 
tions are beginning to be performed, e.g. in fig. [13| we show the lattice spacing obtained by the 
SESAM and T%L collaborations from four physical quantities in both quenched and unquenched 



simulations [36]. In the quenched case there is a spread of results of about ±10%, whereas in the 
unquenched case the spread is smaller (although the errors are still sizeable for some of the quan- 
tities used). In the next 3-5 years it should be possible to compute most of the physical quantities 
discussed below without imposing this approximation. 



2.3.2 Lattice Calculations of B — * p,ir Semileptonic Decays: 

Having discussed the basics of lattice computations, I now briefly present some results from recent 
simulations studying B — > p, 7r semileptonic decays. I start by making the simple observation that from 
lattice simulations we can only obtain the form factors for part of the physical phase space for these 
decays. In order to control discretization errors we require that the three-momenta of the B, ir and p 
mesons be small in lattice units. This implies that we determine the form factors at large values of 
momentum transfer q 2 = (ps — Pn,p) 2 - Experiments can already reconstruct exclusive semileptonic 
b — » u decays (see, for example, the review in []2"ij]) and as techniques improve and new facilities 
begin operation, we can expect to be able to compare the lattice form factor calculations directly with 
experimental data at large q 2 . A proposal in this direction was made by UKQCD p7| ] for B° — > p + l~i>i 
decays. To get some idea of the precision that might be reached, they parametrize the differential decay 
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Fig. 14: Differential decay rate as a function of q 2 for the semileptonic decay B° — > p + l~ui, taken from | p7[ |. Points are 
measured lattice data, solid curve is fit from eq. ( |S5| ) with parameters given in eq. (^). The dashed curves show the variation 
from the statistical errors in the fit parameters. The vertical dotted line marks the charm endpoint. 



rate distribution near g max by: 

dr(B° -> p+l-ut) 12 CP^tf 2 1 2 2 / 2 2 N 

= 1 ° 192vr3Ml g A2({?) ° ( 1 + 6 (ff-w)J. (65) 

where a and 6 are parameters, and the phase-space factor A is given by X(q 2 ) = {m 2 B + m 2 — q 2 ) 2 — 
Am 2 B m 2 . The constant a plays the role of the IW function evaluated at u = 1 for heavy-to-heavy 
transitions, but in this case there is no symmetry to determine its value at leading order in the heavy 



quark effective theory. UKQCD obtain Q37p 

a = 4.6+ g;| ± 0.6 GeV and b = (-8+ 1) x 10~ 2 GeV 2 . (66) 

The fits are less sensitive to b, so it is less well-determined. The result for a incorporates a systematic 
error dominated by the uncertainty ascribed to discretization errors and would lead to an extraction of 
\V u b\ with less than 10% statistical error and about 12% systematic error from the theoretical input. The 
prediction for the dF/dq 2 distribution based on these numbers is presented in Fig. 14. With sufficient 



experimental data an accurate lattice result at a single value of q 2 would be sufficient to fix \ V u b\. 

In principle, a similar analysis could be applied to the decay B° — ► 7r + Z _ P;. However, UKQCD 
find that the difficulty of performing the chiral extrapolation to a realistically light pion from the unphys- 
ical pions used in the simulations makes the results less certain. The B —> 7r decay also has a smaller 
fraction of events at high q 2 , so it will be more difficult experimentally to extract sufficient data in this 
region for a detailed comparison. 

We would also like to know the full q 2 dependence of the form factors, which involves a large 
extrapolation in q 2 from the high values where lattice calculations produce results, down to q 2 = 0. In 
particular the important radiative decay B — ► K*^ (which is dominated by penguin diagrams in the 
Standard Model) occurs at q 2 = 0, so that existing lattice simulations cannot make a direct calculation 
of the necessary form factors. Much effort is being devoted to extrapolating the lattice results down to 



small values of q 2 exploiting the HQET, light-cone sum-rules J38| ] and axiomatic properties of quantum 



field theory. These studies are beyond the scope of these lectures and I refer the students to the review 



articles Q39[ , |40[ ] for a discussion and references to the original articles. 



3 Lecture 3: K° - K° and B° - B° Mixing 

In this lecture I will discuss K-K and B-B mixing in turn, and consider the implications of the experi- 
mental measurements of these processes for our understanding of the unitarity triangle. 



3.1 K°-K° Mixing, e K and e'/e 

The process of K°-K° mixing appears in the Standard Model at one-loop level through the box-diagrams 
of fig. U (with the 6-quark replaced by the s-quark), together with possible long-distance contributions. 
The CP-eigenstates (K\ and K 2 ) are linear combinations of the two strong-interaction eigenstates fj[ 

\K X ) = i= (|^°) + \K )) CP\Kt) = \K\) (67) 

and 

\K 2 ) = -L (\K°) - \K )) CP\K 2 ) = -\K 2 ) . (68) 

Because of the complex phase in the CKM-matrix, the physical states (the mass eigenstates) differ from 
\K\) and \K 2 ) by a small admixture of the other state: 

\K^)+e\K2) \K 2 ) + e\K 1 ) ^ 

^ = , and \ R L> = , 7^^~ ' (69) 

(l + |e| 2 )2 (l + |e| 2 )2 

(the parameter e depends on the phase convention chosen for \K°) and \K )). 

For exclusive decays of iv'-mesons (which have angular momentum zero) into two or three pion 
states, the two pion states are CP-even and the three-pion states are CP-odd. This implies that the 
dominant decays are: 

K$ — > vr7r and Kl — ► 37r , (70) 

which is the reason why Kl is much longer lived than K$ (L and S stand for long and short respectively). 
Kl and K$ are not CP-eigenstates, however, and the decay Kl — > 2tt may occur as a result of the 
small component of \K\) in \K£), and similarly the decay K$ — > 37r may occur because of the small 
component of \K 2 ) in \Kg). Such CP-violating decays which occur due to the fact that the mass 
eigenstates are not CP-eigenstates are called indirect CP -violating decays. As mentioned above, the 
parameter e depends on the convention used to define the phases of the states. A measure of the strength 
of indirect CP-violation is given by the physical parameter defined by the ratio: 

A (Kl — * (im)i— n) , o 
€K = aTtF^ f" = (2.280 ±0.013) 10- 3 e^ , (71) 

where A represents the amplitude for the corresponding process and 1 = implies that the two pions are 



in an isospin state. The numerical value in eq. ( |71[ ) is the empirical result. 



Directly CP-violating decays are those in which a CP-even (-odd) state decays into a CP-odd 
(-even) one: 

K L oc K 2 + IK X . (72) 
Direct (e') 



— ► TTTT 

Indirect (ex) 

TTTT 



Consider the possible quark subprocesses contributing to K — ► -kit decays shown in fig. 15 . The diagrams 



in fig. |l5|(b) and (c) are purely real, whereas that of fig. 15(a) is complex. The two pions in the final state 

8 I use the phase convention so that CP\K°) = \K°) . 



can have total isospin I = or 2; I = 1 is not possible by Bose symmetry. The final states in the 



diagrams of fig. |15|(a) and (b) clearly have 7 = 0, whereas that in the diagram of fig. |15j(c) can have 
either I = or I = 2. Thus direct CP- violation in kaon decays manifests itself as a non-zero relative 
phase between the 1 = and 1 = 2 amplitudes. 




d s 



d d 



u s 



u d 



u 




1 = 0, Complex 
(a) 



1 = 0, Real 
(b) 



I = or 2, Real 
(c) 



Fig. 15: Three subprocesses contributing to if — » tttv decays. 



Of course the gluonic corrections to the diagrams in fig. [l5| generate strong phases, which we de- 
note by 5o,2> where the suffix or 2 represents the isospin of the final state. These phases are independent 
of the form of the weak Hamiltonian. Using the Clebsch-Gordan coefficients for the combination of two 
isospin 1 pions we write: 



A(K° - 
A(K° 



7T 7T 



vrV) 



A e iSo + \I-A 2 e i5 * 



2 . ; A I 

2 ■ iS _ 2 I 1 



A e 



A 2 e ld2 



3 V 3 

The parameter e', which is used as a measure of direct CP-violation in kaon decays, is defined by: 



where [] 



-> = — 

Re A 2 
Re An 



j4 



and 



Im A 2 Im A 



Re A 2 Re A 



- + 5 2 - 6 ~ - 



(73) 
(74) 

(75) 
(76) 



Experimentally the two parameters ex (which, following standard conventions I rename from now 
on as e, e = ejc) and e' can be determined by measuring the ratios: 



V+- = 



A(K L 



7r°vr ) 



A(K S 
A(K L 



7T°7r ) 



7T + 7r 



A(K S -> 7T+7T-) 



2e' 



e + e'. 



(77) 



(78) 



Direct CP- violation is found to be considerable smaller than indirect violation. By measuring the decays 
and using 



??oo 



V+- 



~ 1 — 6 Re 



+ 



(79) 



the NA31 0] and E371 |Kp experiments find (23 ± 7) 10" 4 and (7.4 ± 5.9) respectively for e'/e. 
The Particle Data Group [ 7j] summarises these results as: 



e'/e= (1.5 ±0.8) 10" 



(80) 



e' is manifestly zero if the phases of the 1 = and 7 = 2 amplitudes are the same. 



We will have to wait for the results of the current generation of experiments, which will have a sensitivity 
of O(10~ 4 ), to confirm or deny whether the value is indeed different from zero. 



3.1.1 Determination of e 

In this subsection I discuss the theoretical determination of e, following the discussion and notation in 
the review article by Buras and Fleischer [[I]]. We need to know the matrix element: 



(K°\n^f 2 \K°) 



(81) 



where we write the effective Hamiltonian in the form: 

G 2 



n_iAS=2 



167T 2 



F M 2 



X c 7]iSo{x c ) + X t 7]2S {x t ) + 2X c X t r] 3 So(x c ,xt) 



1 , «i 3 V) ' 



O 



AS=2 



{fx) + h.c. , 



(82) 



the renormalization scale /j is chosen to be much smaller than m c , x\ = mj/Myy and Aj = VidV*. (the 
dependence on X u is eliminated by unitarity, x u = 0). The non-perturbative QCD effects for K°-K° 
mixing are all contained in the matrix elements of the single local composite operator 



O 



AS=2 



(83) 



the remaining terms in eq. ( |82| ) can be calculated in perturbation theory. The So's are the expressions for 
the box-diagrams without QCD corrections: 



S (x) 



4x — lis 2 



4(1 - x) 



S (x c ,x t ) = x c 



In 



3s 3 ln(x) 
" 2(1 -xf 

3x t 3x 2 ln(x t ) 



+ ■■■ ■ 



4(1 -s t ) 4(1 -s t ) 

7/1 , f]2 and r]3 are mass-dependent short-distance QCD factors, at NLO (next-to leading order) 
rn = 1.38 ± 0.20 772 = 0.57 ± 0.01 and 773 = 0.47 ± 0.04 , 



(84) 
(85) 

(86) 



where the errors reflect the uncertainty in Aqcd and the quark masses. The remaining factors in eq. (82) 
are short-distance QCD effects, whose /x-dependence cancels that in the operator O as=2 (in the NDR- 
renormalization scheme J3=1.895). 

As mentioned above, all the non-perturbative QCD corrections are contained in the matrix ele- 
ments of A5=2 . It is conventional to introduce the Bk parameter by the definition: 



{K» I sy*(l - f)d J7„(l - |K U ) = -m^BKitJ,) 



(87) 



The motivation for introducing the parameter Bk in this way comes from the vacuum saturation approx- 
imation (in which Bk = 1), but there is no loss of generality in this definition. Bk{h) depends on the 
renormalization scale, and it is convenient to introduce an (almost) renormalization group invariant Bk 
by: 



B 



K 



1 , jV) ' 

1 + ^r Js 



(88) 



The weakest link in the theoretical calculation of e is the evaluation of Bk- Many lattice computations 
have been performed, and a recent compilation of results gave Bk = 0.90 ± 0.06 [43]. A calculation 



based on the 1/N approximation gave Bk = 0.70 ± 0.10 [44]. Buras and Fleischer perform their 
analysis with 

B K = 0.75 ±0.15 . (89) 
e is given in terms of the quantities introduced above by 

e = C e B K Im \ t {Re A c [?7i5o(x c ) - rj 3 So( x c a:*)] - Re A t r&Sofct)} e'* , (90) 



where 



c = FJK W = 3 J8 . 1Q 4 

6V2tt 2 (AM k ) 



(91) 



and the mass difference AMk = Mr l — Mk s - 

We obtain information about the unitarity triangle by comparing the theoretical prediction with 
the experimentally measured value of e. In terms of fj and p the prediction is the hyperbola 



V 



(1 - p)A 2 m S (x t ) + P 



A 2 B K 



0.226 



(92) 



where P = 0.31 ± 0.02. This is illustrated in fig. |T§ 




C(0,0) 



B(l,0) 



Fig. 16: A precise determination of e would fix the vertex A to lie on a a hyperbola (schematically represented by the solid 
curve). 



3.1.2 The AI =1/2 Rule and e'/e 

In this subsection I briefly review the elements in the theoretical predictions for e'/e and related quanti- 
ties. I start by rewriting eq. (^) in the form: 

<? = ~Z{l-n)e i *, (93) 

where 

ReAo ImAo „ 1 Im^o tt it 

W = 1T^< f = 5-T> Q = -T^ and 0=- + 5 2 -5 o ^-. (94) 
Re Aq RcAq ulmAo 2 4 

Even after about 40 years, we still do not understand theoretically the AI = 1/2 rule, i.e. the 
empirical observation that u ~ 1/22. The rate for transitions in which the isospin changes by 1/2 
(AI = 1/2 transitions) is enhanced by a factor of about 450 over those in which it changes by 3/2 
(AI = 3/2 transitions). There is a relatively small enhancement of the AI = 1/2 amplitudes, by 
a factor of 2-3, due to perturbative effects (i.e. in the calculation of the Wilson coefficient function). 
The remaining enhancement is assumed to be due to long-distance QCD effects, perhaps because the 
matrix-elements of penguin operators are large, but this has still to be convincingly demonstrated. 



In the recent analyses of e'/e, one combines the experimental values for the real parts of the 
amplitudes (Re A = 3.33- 1CT 7 GeV, Re A 2 = 1.50- 1(T 8 GeV, u = 0.045), with theoretical predictions 
for the imaginary parts. The effective Hamiltonian for these A5 = 1 transitions takes the form: 

(-t 10 

n eff (As = i) = -| Yl + T v^)) . (95) 

where r = —VtdV t */V u dV* s . There are 10 independent operators in the AS = 1 effective Hamiltonian, 
and the major uncertainty is in the values of their matrix elements and in the mass of the strange quark 



m s . The most complete analyses have been performed by Buras et al [45] who find: 

(-1.2 • 10~ 4 ) < e'/e < 16 • 10~ 4 (/) ; e'/e = (3.6 ± 3.4) 10~ 4 (II) , (96) 
where / and II refer to two different ways of treating the many uncertainties (more and less conservative 



respectively), and by the Rome group [ |46| ] who find 

e '/e = (4.6 ± 3.0) 10" 4 . (97) 

The various contributions to e'/e have different signs, and there are significant cancellations. For this 
reason even the expected sign cannot be predicted with full confidence. 

I end this subsection with a brief summary of the main points: 

• An understanding of the AI = 1/2 rule would be an important milestone in controlling non- 
perturbative QCD effects. This is a realistic, but non-trivial, challenge for the lattice community. 

• A measurement of a non-zero value for e'/e would be a very important qualitative step in particle 
physics; confirming the existence of direct CP-violation. 

• The uncertainties in the values of the matrix elements, make it difficult to make a precise prediction 
for e'/e. We expect it to be at the several xlO -4 level, but accidental cancellations between 
contributions with opposite signs may make it smaller than this. In the coming years, experiments 
at CERN, FNAL and at Da^ne will be sensitive to a value of about 10~ 4 . 

3.2 B°-B° Mixing 

In this subsection we consider B-B mixing, from which we get information about the V t d and Vt s ele- 
ments of the CKM-matrix. I start with some general formalism for the mixing of a neutral psudoscalar 
meson P° with its antiparticle P°. I write the wave-function in two-component form: 



\m)= (Jg I ^a(t)\P»)+b(t)\P»). (98) 
The time dependence of \ip(t)) is given by the Schrodinger equation: 

if t \m) = ( M - y ) IVW> , (99) 

where M — iT/2 is the 2 x 2 "mass-matrix", whose elements are equal to (P®\H e ff\Pj) /2Mp. Using 
the optical theorem, the absorbtive part of this matrix comes from real intermediate states: 

T a = 2M; 5D< p i° \ Hw |n)(n \ Hw |jP > 0) 27r 5{ - En ~ Mp) ■ (100) 



For B°-B° mixing, there are no contributions to Ti2 from intermediate states containing the top quark. 



The mass-matrix takes the form 



M 



T 




(101) 



where the diagonal terms are equal by CPT-invariance, and A, p 2 and q 2 are complex. The eigenvalues 
of this matrix are Mi 2 = A ± pq, so that the mass difference between the two physical eigenstates is 



AM = Mi - M 2 = 2pq = 2(Mi2 - iTi 2 )^(M^ 2 - iT* 12 ) 



(102) 



For the i?-system, the box diagrams have contributions from t-quarks in the loops (see fig. |J). 
By unitarity these contribute only to the real part of the mass-matrix. The mass behaviour of the box 
diagrams in fig. and the large value of V t b imply that these contributions dominate and that | Ti2 1 <C 
IM12I for B-B mixing (this is not the case for K-K mixing). Thus AM ~ 2|Mi 2 |, and can be calculated 
from the box diagrams. 

The evaluation of the box diagrams follows similarly to that in the kaon system. The non- 
perturbative QCD effects are contained in the matrix element of the AB = 2 operator: 



O 



AB=2 



67^(1 -7 5 )d67 M (l- 7 5 )d 



(103) 



and it is also convenient and conventional to introduce the Bb parameter analogously to the definition of 



Bk in eq. (|87[). The theoretical expression for the mass difference is: 



AM 



VB M B B b f 2 B M 2 , S (x t ) \V td \ 2 



(104) 



where tjb is a pertubative QCD correction (tjb = 0.55 ± 0.01). Throughout this discussion I have been 
implicitly assuming that we are considering neutral B^ mesons, but an analogous discussion holds for 
the B s system. In the kaon system, all the non-perturbative QCD effects were contained in the parameter 
Bk, since the leptonic decay constant fx is a measured quantity. For the i?-mesons this is not the 
case, and one is obliged to use model estimates or lattice or sum-rule calculations for both the decay 
constants fs and the Bb -parameters (or of the matrix element of AB=2 itself). For example, in a recent 
compilation of lattice results we obtained [[^] : 



f B . = (170 ± 35) GeV, f Ba = (195 ± 35) GeV, B Bd = 1.4(1) 



and 



Jb s \J B Bs 



fB d ^JBB d 
Combining all the elements one obtains: 



1.14(8) . 



AM Bd = 0.50 ps" 1 
AM Bs = 15.1 ps" 1 



200 MeV 



1.52 



\V t d\ rj B 



240 MeV 



ni t (m t ) \ 
170 GeV/ 0.0088 0.55 



/ m t (m t ) \ L52 \V ts \ rj B 
V 170 GeV J 0.040 0.55 ' 



(105) 
(106) 

(107) 
(108) 



The mass difference AM is a measure of the oscillation frequency to change from a B° to a B° 
and vice-versa. Imagine that we start with a B° at time t = 0, then the propabilities that we have a B° 
or B° at a later time t are exp(— Tt) cos 2 (AMt/2) and exp(— Ft) sin 2 (AM t/2) respectively. Defining 
P{B ) and P(B°) to be the probabilities (obtained by integrating over time) that when the meson decays 
it is a B or B respectively, we obtain 



P(B C 



+ 



r r 2 + (AM) 2 



and P(B C 



r 2 + (am) 2 



(109) 



The ratio of these integrated probabilities is given by 

_ P(B°) x 2 

where x = AM/T. The 1996 Particle Data Group review quotes: 

AM Bd = (0.474 ± 0.031) ps" 1 and x d = 0.73 ± 0.05 
in agreement with theoretical expectations and 

AM Bs > 5.9ps _1 and x s > 9.5 at 95% CL . 



(110) 



(111) 



(112) 



The oscillation in the B s system are too rapid for observation in experiments which have been performed 
up to date. 

What do these measurements imply for the determination of the position of the unitarity triangle? 
In principle a measurement of AMs d allows for a determination of the Vtd element of the CKM matrix. 
The side BA of the unitarity triangle is given by: 



\v td v t i\ _ 1 \v td \ 



(1 - p) 2 + f) 2 



(113) 



\V cd V* b \ A \Vcb\ 

so that a measurement of V t d (assuming knowledge of V&) would imply that the locus of possible posi- 



tions of the vertex A lies on a given circle centred on B, see fig. 17 




C(0,0) 



5(1,0) 



Fig. 17: A precise determination of \V t d/V c b\ from studies of B — B mixing would fix the vertex A to lie on a circle centered 
on the vertex B (schematically represented by the solid curve). 



4 Lecture 4: CP- Violation in B-Decays and Inclusive U-Decays 

In this final lecture I will consider two topics in B -physics, CP- violation (and mixing-induced CP- 
violation, in particular) and the application of the heavy quark expansion to inclusive decays (and to 
beauty-lifetimes and the semileptonic branching ratio, in particular). I will also make some brief remarks 
about how difficult it is to make predictions for exclusive nonleptonic decays. 

4.1 CP- Violation in B -Decays 

Measurements of CP-asymmetries in neutral P-meson decays into CP-eigenstates will allow us to 
address 3 fundamental questions: 



i) Is the phase in the CKM matrix the only source of CP-violation? 



ii) What are the precise values of the CKM parameters? We shall see that the principal systematic 
errors, particularly those due to hadronic uncertainties, are reduced significantly in some cases. 



iii) Is there new physics in the quark-sector? 



A particularly important class of decays of neutral P-meson systems, are mixing induced CP- 
violating decays (these are not possible in decays of charged P-mesons). I now briefly review this 
important topic, for which we need to return to the formalism introduced in the discussion of P-P 



mixing in subsec. above. The two neutral mass-eigenstates are given by 

1 / ,„n, . ,„ . 1 



\Bl) 



VP 1 



p|P u ) + g|P u ) and \B H ) 



\Jp 2 + q 2 



p\B u ) -q\B u 



(114) 



Starting with a B° meson at time t = 0, its subsequent evolution is governed by the Schrodinger equation 



( 



\B° phys (t)) = 9 + (t)\B°) + (?) g-(t)\E°) , 



where 

9+{t) = ex P 
flL-(i) = exp 
and M = (M H + M L )/2 0. 



Tt 
Y 
Tt 

T 



exp[— iMt] cos 
exp[— iMt] i sin 



(AMt 
\ 2 
/AMt 
\ 2 



(115) 

(116) 
(117) 



4.1.1 Decays of Neutral B -Mesons into CP-Eigenstates 
Let fcp be a CP-eigenstate and A, Abe the amplitudes 

A = (fcp\n\B°) and A = (f C p\H\&°) 

Defining 



A 



q A 
p A 



(118) 



(119) 



we have 



(f C P | H | P^ hys ) = A [ 9+ (t) + A 5 _(t)] and (f CP \ H \ B» hys ) = A V - \g_(t) + Xg+(t)] . (120) 

The time-dependent rates for initially pure B° or B° states to decay into the CP-eigenstate fcp at time 
t are given by: 



T(B° phys (t) ^ fcp) = \A\ 2 e 



2 -Ft 



T(B°(t)^fcp) = \A\ 2 e- rt 



1 + J A ^ 2 + 1 2 * A * 2 cos(AMt) -ImA sin(AMt) 
1-|A| 2 , 1-|A| 2 



+ 



cos(AMt) +ImA sin(AMt) 



The time-dependent asymmetry is defined as: 

r(P p ° hy sW - fcp) - r(B° phys (t) fcp) 



r(£° hys (t) - fcp) + r(B% ys (t) f C p) 

(1 - |A| 2 )cos(AMt) - 21mA sin(AMt) 
1 + IAI 2 ' 



(121) 
. (122) 

(123) 
(124) 



"Starting with a B° meson at t = 0, the time evolution is \B° hys (t)) = (p/q) g-(t)\B°) + g+(t) \B°). 



If, as is the case in some important applications, \q/p\ = 1 (which is the case if Ar << AM) and 
\A/A\ = 1 (examples of this will be presented below), then |A| = 1 and the first term on the right-hand 



side of eq.(|124j) vanishes. 

The form of the amplitudes A and A is: 

A = A i eiSl e * and A = J2 A i eiSl ( 125 ) 

i i 

where the sum is over all the contributions to the process, the A4 are real, Si are the strong phases and 
the <f>i are the phases from the CKM matrix. In the most favourable situation, all the contributions have 
a single CKM phase (<j>£> say) and A/ A = exp(— 2i4>d)- Under the assumption we are making that 



T 12 « M12, q/p = ^M{ 2 /M l2 = exp(-2i</> M ), and A = exp(-2i(0D + <f> M )). Thus 

ImA = -sm(2(0 D + M )). (126) 
From the box diagrams of fig. |5| we obtain: 

^ -WS and (*) = (127) 



p;B d y t *d v tb \pJb s v&vtb 

To illustrate the above discussion let us consider three processes in which the 6-quark decays 
through the subprocess b — > djUiUi. The corresponding tree-level diagram is 




(128) 



B d — ► J/* K s : In this case 

\(B -> J/^K S ) 



V td V* V cs V* d V cb V c 



V* d V tb V*V rA V1V, 



cd V cb Vcs 



sin(2/3) 



(129) 



The first factor in eq. (129) is the factor (q/p) B d m ec l- (127). the third factor is the analogous one for the 
final state kaon, and the second factor is A/A as in eq. (128) with m = c and bj = s. We recall from the 
discussion of the first lecture that the angle 



j3 = arg 



V c dV. 



cb 



\ V td V* b ) 

There is also a small penguin contribution to this process, due to the subprocess: 



(130) 




which is proportional to Vt b V t *, whose phase is, to a very good approximation, equal to that of V cb V* s . 
Thus hadronic uncertainties are negligible in the determination of the angle from this process, and for 
this reason we can consider it a gold-plated one. This is an (almost) ideal situation. 



Bd — > 7r + 7r : For this process, in the notation of eq. ( 128), we have U{ = u and dj = d, so that 



a v ub v: d { } 

and Im A = sin(2a), where from the first lecture we recall that 

-MM)- 

In this case there is also a small penguin contribution, which is proportional to V td V t b, which has 
a different phase from the tree diagram. This implies that the hadronic uncertainites are greater than for 
the process B d — > J/^Ks, although they are still reasonably small (and can be reduced further using 
isospin analysis). 



B s — ► p °Ks : For this process in the notation of eq. (128), we have m = u and dj = d again, so that 

v ts v t t v ub v: d v cd v* 



\(B S -> pK s ) 



Using the unitarity relations and the fact that 



v t * s v tb v: b v ud v c * d v cs 



(133) 



7 = arg 



(134) 



we find that Im A = — sin(27). 

Although this process is frequently proposed as a potential way of measuring the angle 7, the 
penguin contributions are relatively large and Buras and Fleischer stress that it is a wrong way [[]]]. It has 
been suggested that it might be possible to use combinations of B — > DK amplitudes to determine the 
angle 7 p7[ ] . This is illustrated in fig. 18, where | D+ ) = 1/ v2( | D ) + 1 D°) ) , but clearly a measurement 
of all the amplitudes is a huge experimental challenge. 




Fig. 18: Six amplitudes from which the angle 7 can (in principle at least) be determined. 



The exciting experimental and theoretical program of research into mixing induced CP-asymetries 
will continue at the (near) future B-factories and dedicated ^-experiments. 



4.1.2 CP-Violation in Charged B-Decays 

Observation of CP- violation in decays of charged P-mesons would demonstrate the existence of direct 
CP-violation, but the hadronic uncertainties make it unlikely that they can be used to determine the 
parameters of the unitarity triangle with the required precision. Since we need an interference between 
different contributions to the amplitude, we write the amplitudes in the form: 

A(B~ -> f) = V x A 1 e iSl + V 2 A 2 e i&2 and A(B + -> /) = V{ A x e iSl + V 2 * A 2 e l52 , (135) 

where I have exhibited two contributions to each amplitude, V\ j2 are CKM-matrix elements and the 
Si j2 are strong interaction phases. For example in the process B + —>■ ir°K + the interference is between 
curent-current and penguin contributions, and in the decay B + — ► K°K + it is between penguin diagrams 
with internal it and c quarks. 

The CP- violating asymmetry is now defined by 

_ T(P+ -> f)-T(B- f) _ 21m (Vjjg) sin(fr-£ 2 ) A X A 2 

^ CP r( B + _ /) + F (p- _ /) \V 1 \ 2 A 2 + \V 2 \ 2 A 2 + 2R e (V 1 V 2 *)cos(S 1 -S 2 )A 1 A 2 ' 1 ' 

from which we see that in order to have a non-zero asymmetry we require both Im (VxV 2 ) and sin(<5i — 
6 2 ) to be non-zero. The strong interaction phases are very difficult to quantify, so it is hard to use 
measurements of these CP-asymetries to determine the CKM matrix elements. 

4.2 Inclusive Nonleptonic Decays of Beauty Hadrons 

In this subsection I discuss inclusive nonleptonic decays of beauty hadrons in general, and two very 
interesting phenomenological problems in particular the lifetimes of the hadrons and their semileptonic 



branching ratios. The discussion will use the formalism of Bigi et al.(see ref. Q48p and references therein), 
developed and used by them and many other groups, in which inclusive quantities are expanded in inverse 
powers of the mass of the heavy quark, e.g. 



r(i?b) = \c 3 [l + ^P^ + c 5 ^| + 0(-A,U , (137) 




G 2 F ml\V ch \ 2 J / A! + 3 A 2 ^ , A 2 
192tt3 \ C3 { 2m 2 

where F is the full or partial width of a beauty hadron H^, C3 5 are coefficients which can be computed 
in perturbation theory (the coefficient functions obtained when matching QCD onto the HQET) and Ai^ 
are the matrix elements of the kinetic energy and chromomagnetic operators respectively: 

X 1 (B) = ^-(H b \mD) 2 h\H b ) and A 2 = \ (J^H^W {m 
2m Hh 3 2mm 



and h is the field of the heavy (static) quark. An important feature of the general expression in eq. (137) is 
the absence of terms of 0(l/mb), which is a consequence of the absence of any operators of dimension 
4 which can appear in the corresponding OPE [0^] . 



I will not discuss the derivation of the expansion in eq. ( |137[ ) in detail. One starts, similarly to the 
derivation of the OPE for moments of deep-inelastic structure functions, by using the optical theorem 
which states that the width is given by the imaginary part of the forward elastic amplitude. In this 
case this amplitude has two insertions of the weak Hamiltonian. Since the 6-quark is heavy and decays 



into lighter states, the separation of these two insertions is small, and eq. ( |137| ) is the corresponding short 
distance expansion (for the structure functions it is a light-cone expansion). We now apply this expansion 
to a study of beauty lifetimes and the semileptonic brabching ratio. 



4.2.1 Beauty Lifetimes 



Using the expression in eq. ( |137| ) for the widths, one readily finds the following results for the ratios of 
lifetimes: 



r(B-) 
t(B°) 

r(A 6 ) 
t(B°) 



l + O 



1 



mr 



I + *4(A») -ffii?) + ^ /4(A fe )-^(ff) + Q (l^ 



2m\ 



mi 



(0.98 ±0.01) + [ —5 



(139) 



(140) 



where = — and [i G = 3A2- In order to obtain the result in eq.(|140J), one needs to know the 
difference of the kinetic energies of the 6-quark in the baryon and meson. To leading order in the heavy 
quark expansion we have: 



/4(Afc)-A4(B) 



M B M D ( M Ab - M Ac 
2 V M B - M D 



3 M B * - M D * 

4 M B - M d 



(141) 



From equation (141), and using the recent measurement of m^ B from CDF [|50|], one finds that the 
right hand side is very small (less than about 0.01 GeV 2 ). The matrix elements of the chromomagnetic 
operator are obtained from the mass difference of the B*- and S-mesons (see eq.( 138[ )) and from the 



fact that the two valence quarks in the A& are in a spin-zero state (which implies that //g,(A&) = 0). The 



theoretical predictions in eqs. (|13 

r(B-) 
t(B°) 



and (|140J) can be compared with the experimental measurements 



1.06 ±0.04 and 



t{B°) 



0.79 ± 0.05 



(142) 



The discrepancy between the theoretical and experimental results for the ratio r(Af,)/r(i3 ) in eqs. (140) 
and ( |142j ) is notable. It raises the question of whether the 0(1/ ml) contributions are surprisingly large, 
or whether there is a more fundamental problem. I postpone consideration of the latter possibility and 
start with a discussion of the 0(l/-m|) terms. 
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Fig. 19: Examples of diagrams whose imaginary parts contribute to the total rates for the decays of beauty hadrons (left-hand 
sides) and the operators they correspond to in the Operator Product Expansions. V represents a Dirac matrix. 



At first sight it seems strange to consider the 1 / corrections to be a potential source of large 
corrections, when the 0(l/m|) terms are only about 2%. However, it is only at this order that the "spec- 
tator" quark contributes, and so these contributions lead directly to differences in lifetimes for hadrons 
with different light quark constituents (consider for example the lower diagram in Fig. 19, for which, 



using the short-distance expansion, one obtains operators of dimension 6). Moreover, the coefficient 
functions of these operators are relatively large, which may be attributed to the fact that the lower dia- 
gram in Fig. [l9| is a one-loop graph, whereas the corresponding diagrams for the leading contributions 
are two-loop graphs (see, for example, the upper diagram of Fig. [l9|). The corresponding phase-space en- 
hancement factor is \6ir 2 or so. We will therefore only consider the contributions from the corresponding 
four-quark operators, neglecting other 0(l/m%) corrections which do not have the phase space enhance- 
ment [MJ]. For each light-quark flavour q, there are four of these |^]: 

Ot =6 7m (1 - 7 5 )g q^(l -7 5 )& ; 2 = 6(l+ 7 5 )gg(l + 7 5 )6 (143) 
T 1 = 6 7M (l-7 5 )r a g g 7 "(l - 7 5 )T a 6 ; T 2 = 5(1 + j 5 )T a q q{l + 7 5 )T a 6 (144) 

where T a are the generators of colour SU(3) and q represents the fields of the light quarks. Thus we 
need to evaluate the matrix elements of these four operators. 



For mesons, following ref. [51J, I introduce the parametrization 

{B\O t \B)\^ =mb = BiflMl ■ {B\Ti\B)\^ =mb = e t f 2 M 2 , (145) 
where /i is the renormalization scale. We have chosen to use as the renormalization scale. Bigi 



et al. (see ref. [ |48| ] and references therein) prefer to use a typical hadronic scale, and estimate the 
matrix elements using a factorization hypothesis at this low scale. Operators renormalized at different 
scales can be related using renormalization group equations in the HQET (sometimes called hybrid 
renormalization [|2|]). For example, if we assume that factorisation holds at a low scale \x such that 
a s (/j, 2 ) = 1/2, then, using the (leading order) renormalization group equations, one finds B\ = B 2 = 
1.01 and ei = €2 = —0.05 0. In the limit of a large number of colours N c , Bi = 0(N®) whereas 
e f = 0(l/JV c ). 

For the heavy quark symmetry implies that 

(A 6 |0 2 |A 6 ) = -i(A b |0 1 |A 6 ) and (A 6 |T 2 |A 6 ) = -i(A 6 |:Zi|A 6 ) , (146) 

so that there are only two independent parameters. It is convenient to replace the operator T\, by 0\ 
defined by 

Ox = 6*7m(1 " 7 5 V - 1 5 W , (147) 

where i, j are colour labels, and to express physical quantities in terms of the two parameters B and r 
defined by 

(AftlOilAb)^ = -S(A 6 |Oi|A & ) M=rnb (148) 
' -(AalOilAj)^ = - f -^-r. (149) 



2M Ab x U| X| - 48 

We do not know the values of these parameters. In quark models B = 1, and r = 0.2-0.5. Using 
experimental values of the hyperfine splittings and quark models, it has been suggested that r may be 
larger [||], e.g. 

^It^ = ^±0.1 and r ^ = 1.8 ± 0.5 . (150) 



"i use the notation of ref. 
12 By factorizatii 
and respectively. 



12 By factorization we mean that if the Bi's and e;'s had been defined at this scale (instead of rnb) they would have been 1 



The lifetime ratios can now be written in terms of the six parameters B12, e\^, B and r (as well 

as /b): 



t{B- 



t(B°) 
t(B°) 



1 + 



f B 



0.98 + 



200 MeV 
fB 



{0.02Si + O.OOB2 - 0.70ei + 0.20e 2 } 
j {-O.OOSa + 0.00^2 - 0.17ei + 0.20e 2 



200 MeV 
+ (-0.01 - 0.02S) r} , 



(151) 



(152) 



where the effective weak Lagrangian has been renormalized at n = m&. The central question is whether 
it is possible, with "reasonable" values of the parameters, to obtain agreement with the experimental 



numbers in eq. ( |142| ). At this stage in our knowledge, the answer depends somewhat on what is meant by 
reasonable. For example, Neubert, guided by the arguments outlined above, has considered these ratios 



by varying the parameters in the following ranges 054Q : 



B h B G 



2 4 

3' 3 



e* G 



1 1 

3' 3 



r G [0.25,2.5] ; 



fB \ 



200 MeV 



J G [0.8,1.2] . (153) 



He concludes that, within these ranges, it is just possible to obtain agreement at the two standard deviation 
level for large values of r (r > 1.2) and negative values of e 2 . Lattice studies of the corresponding matrix 
elements are underway; a recent QCD sum-rule calculation has found a small value of r, r ~ 0.1- 



0.3 [55] 



If the lattice calculations confirm that the parameter r is small, or find that the other parameters are 
not in the appropriate ranges, then we have a breakdown of our understanding. If no explanation can be 
found within the standard formulation, then we will be forced to take seriously the possible breakdown 



of local duality. This is beginning to be studied in toy field theories [56 



4.2.2 The "Baffling" Semileptonic Branching Ratio 



This was the name given by Blok et al. [ p8| ] to the observation that the experimental value of the semilep- 
tonic branching ratio 

B SL = T{B - ^ (154) 

Hi r(-B — * xiv) + Thad + r mre 



appeared to be lower than expected theoretically. In eq. ([154]) the sum is over the three species of lepton, 
and T^ad and T rare are the widths of the hadronic and rare decays respectively. Bigi et al. concluded that 
a branching ratio of less than 12.5% cannot be accommodated by theory [|5"8|]. Since then Bagan et al. 
have completed the calculation of the 0(a s ) corrections, and in particular of the b — > ccs component of 
the decay (including the effects of the mass of the charm quark) [^] ; these have the effect of decreasing 
Bsl- With M. Neubert, we used this input to reevaluate the branching ratio and charm counting (n c , the 
average number of charmed particles per i3-decay) [|T|] finding, e.g. 

B SL = 12.0 ± 1.0% (p = m b ) n c = 1.20 + 0.06 (p = m b ) 

B S l = 10.9 ±1.0% {n = m h /2) n c = 1.21 + 0.06 (p = m h /2) . 1 ' 

[i is the renormalization scale and the dependence on this scale is a reflection of our ignorance of higher 
order perturbative corrections. The experimental situation is somewhat confused, see Fig. In his 
compilation at the ICHEP conference last year []60|], Richman found that the semileptonic branching 
ratio obtained from i?-mesons from the T(45) is 

B SL (B) = (10.23 ± 0.39)% , (156) 



13 Note that the rapporteur at the 1997 EPS conference argued that the branching ratio had been overestimated by the LEP 
collaboration [EH]]. 




Fig. 20: Theoretical Prediction (shaded region) of the semileptonic branching ratio and charm counting. The data points are the 
experimental results from high-energy (LEP) and low energy (i.e. at the T(4S) from CLEO) experiments. 



whereas that from LEP is: 

B SL (b) = (10.95 ± 0.32)% . (157) 

The label b for the LEP measurement indicates that the decays from beauty hadrons other than the 
S-meson are included. Using the measured fractions of the different hadrons and their lifetimes, and 
assuming that the semileptonic widths of all the beauty hadrons are the same, one finds: 

B SL (b) = (10.95 ± 0.32)% B SL (B) = (11.23 ± 0.34)% , (158) 

amplifying the discrepancy. It is very difficult to understand such a discrepancy theoretically, since the 
theoretical calculation only involves Fsl (and not F^ad f° r which the uncertainties are much larger). 
In view of the experimental discrepancy, I consider the problem of the lifetime ratio t(Ai ) )/t(B°), 
described in subsection [1.2. 1 above, to be the more significant theoretical one. 



4. 2. 3 Exclusive Decays of B -Mesons 

A large amount of experimental data is becoming available, particularly from the CLEO collaboration 



(see ref. [62] and references therein) on two-body (exclusive) nonleptonic decays of 5-mesons. This is 
an exciting new field of investigation, which will undoubtedly teach us much about subtle aspects of the 
Standard Model. Unfortunately, at our present level of understanding we are not able to compute the 
amplitudes from first principles, and are forced to make assumptions about the non-perturbative QCD 
effects. Frequently these assumptions concern factorization, i.e. whether matrix elements of operators in 
the effective Hamiltonian which are products of two V— A currents can be written in terms of products 
of matrix elements of the currents, e.g. 

(D+7T-\ {du) L (cb) L \ B°) = (D+\ (cb) L \ B°) (tt-| {du) L | 0) , (159) 



where the first matrix element on the right hand side could be obtained from semileptonic decays (see 
sec. 2A) and the second from the known value of f w . In some cases, the factorization hypothesis can be 




motivated by the concept of colour transparency, that a light colour singlet meson, produced with large 
energy at the weak decay vertex (and which is hence a small colour-dipole), will travel a long way from 
the interaction region before it grows to a size where it could interact strongly 54]. It is possible 
that in some cases factorization will give a reasonable first approximation and in many others it will 
not be useful at all. Thus current analyses are limited to a semi-quantitative level, and at this stage it 
is not possible to endorse any of the approaches. I will not discuss these analyses in these lectures, but 
to wish to point out the importance of this emerging field, (see refs. [65] and [ |40| ] for more extensive 
discussions). 



5 Conclusions 



In these lectures I have briefly reviewed the formalism required for the study of the weak decays in the 
Standard Model of Particle Physics. From these studies we are attempting to measure the parameters 
of the Standard Model, to understand the origin of CP-violation and to look for signatures of physics 
beyond the Standard Model. I have explained that the major theoretical difficulty in interpreting the 
wealth of experimental data is our inability to control the non-perturbative strong interaction effects to 
sufficient precision. Througout these lectures I have illustrated the general discussion with applications 
to various physical processes. The current status of the uncertainties in determining the vertex A of the 



unitarity triangle, from the analysis of Buras and Fleischer is reproduced in fig. |21]. The circular arcs are 
obtained from information on B — B mixing, using a modified analysis to that discussed in section 
above. In order to reduce the hadronic uncertainties, we can take the ratio of the expressions for AM^ 
and AMb s in eqs. (107) and (108), and, using the upper bounds on AM^ (0.482 ps _1 ) and |Vf S /V^,| 



(0.993), obtain a bound for the quantity in eq. ( |113[ ): 



l W, 



td\ 



A \V cb \ 



< 



'10.2ps~ 
AMr 



(160) 



where £ is defined in eq. ( 106), where the lattice result is also given. The curves in fig. 21 represent these 



bounds for several choices of £ and AM b s Q 
14 In fig. |l] AM Bb is written as AM S . 



I end these lectures by summarising the tourist guide of Buras and Fleischer ^ to different pro- 
cesses, in which stars are awarded to the processes depending on the level of theoretical uncertainty. 



*** These are processes in which there are no theoretical uncertainties i.e. they are particularly small 
(less than 2% say). Among these is the mixing-induced CP-violating asymmetry in Bd — ► 
J/^Ks discussed in section |4.1.1| , from which we expect to determine the angle of the unitarity 
triangle. Other quantities in this category which I have not had time to discuss include the ratio 
of (inclusive) branching ratios Br(P — > XdfP) /Br(B — > X s vv) from which one would obtain 
l^tdl/l^tsl an d the decays Kl — ► t^vv (which would give Im X t ) and K + — > -k+vv (which would 



give \Vtd\)- For the kaon decays I have explained in section |l.2.l| , that the small mass difference 



between the quarks in the K and tt mesons implies that the hadronic uncertainties are small. 
The processes in this category have small hadronic uncertainties, of the order of 5-10% say. They 



include the evaluation of V c b discussed in sec. 2.1, AM^/ AM S discussed in section 3.2 (from 



which we obtain |Vtd|/|Vt s |) and the mixing-induced CP-asymmetry in Bd — > 7r + 7r~ (from which 
we would obtain the angle a) discussed in section |4.1.1[ There is a large set of other mixing- 
induced asymetries in this category, yielding, in principle, all three angles. Being perhaps a little 
optimistic, the determination of |V^|/|VJ.&| in inclusive decays, also belongs to this category. 

The processes in this category have smallish hadronic uncertainties, perhaps in the region of 15% 
or so. Of the quantities discussed in this course, AMd, AM S (section 3_2) and e (section |3.1.1[ ) 
belong in this category. 

As we have seen in these lectures, there are also important quantities for which the theoretical un- 
certainties are too large for the experimental data to give accurate information about the parameters 
of the Standard Model (and hence these quantities are not ascribed any stars), but which neverthe- 



less give interesting qualitative information. These include e'/e (discussed in section 3.1.2) and 
most CP-asymmetries in decays of charged P's. Non-leptonic decays of P-mesons in general 
belong to this category. 



Of course it is to be hoped and expected that theoretical progress will continue, and that the 
interesting physical quantities will gradually acquire more stars. The theoretical community is working 
very hard to achieve this. 
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